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Abstract 

The rheological properties of fibrin networks have been of long-standing 
interest. As such there is a wealth of studies of their shear and tensile 
responses, but their compressive behavior remains unexplored. Here, by 
characterization of the network structure with synchronous measurement 
of the fibrin storage and loss moduli at increasing degrees of compression, 
we show that the compressive behavior of fibrin networks is similar to that 
of cellular solids. A non-linear stress-strain response of hbrin consists of 
three regimes: 1) an initial linear regime, in which most hbers are straight, 
2) a plateau regime, in which more and more fibers buckle and collapse, 
and 3) a markedly non-linear regime, in which network densification oc¬ 
curs by bending of buckled fibers and inter-fiber contacts. Importantly, 
the spatially non-uniform network deformation included formation of a 
moving “compression front” along the axis of strain, which segregated the 
fibrin network into compartments with different fiber densities and struc¬ 
ture. The Young’s modulus of the linear phase depends quadratically on 
the fibrin volume fraction while that in the densified phase depends cubi- 
cally on it. The viscoelastic plateau regime corresponds to a mixture of 
these two phases in which the fractions of the two phases change during 
compression. We model this regime using a continuum theory of phase 
transitions and analytically predict the storage and loss moduli which are 
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in good agreement with the experimental data. Our work shows that hb- 
rin networks are a member of a broad class of natural cellular materials 
which includes cancellous bone, wood and cork. 


1 Introduction 

Fibrin is a polymer made from the blood plasma precursor protein fibrinogen. 
Fibrin forms a porous network of branching fibers that provides the scaffold of 
protective hemostatic blood clots and pathological obstructive thrombi. The 
structure and properties of the fibrin, including mechanical resilience to forces of 
blood flow and clot contraction, determine the course of pathological conditions, 
such as coronary heart disease, ischemic stroke, or bleeding [1, 2]. Hence, the 
tensile and shear behaviors of fibrin networks have been studied for a long time 
[3, 1]. Recent experimental studies have also focused on the tensile properties 
of a single fibrin fiber including the molecular unfolding of fibrin molecules at 
large deformations [4, 5, 6, 7, 8]. The tensile behavior of single hbrin hbers can 
be used as an input to network models to predict the tensile behavior of the 
networks. We did this recently using the eight-chain model and were able to 
predict the mechanics of hbrin networks under tension [9, 10, 14, 15, 16]. The 
shear behavior of hbrin networks has also been studied both experimentally and 
theoretically [13, 14, 15, 16, 17]. For networks with thin hbrin hbers the shear 
behavior can also be explained in terms of the entropic elasticity of single hbrin 
hbers [14, 16]. 

Even though tension and shear behaviors of hbrin networks have been stud¬ 
ied in detail, compression of these networks remains largely unexplored. Re¬ 
cently, Kim et. al. [18] found that the non-linear stress-strain curves of hbrin 
networks in compression have three distinct regimes, similar to the curves ob¬ 
served in networks of carbon nanotubes [19]. There is an initial linear response 
which is expected of any elastic solid for small strains, followed by a stress 
plateau that begins when hbers in the network start to buckle and continues 
until a majority of them have done so. Then, network densihcation occurs 
creating a large number of contacts between the hbers as the stress increases 
steeply. This tri-phasic mechanical behavior could presumably be modeled us¬ 
ing the theory of foams and other cellular solids [20] as shown schematically 
in Figure 1. The initial shear and compression moduli of foams in the linear 
regime can be computed in terms of the Young’s modulus of the foam material 
and the geometric parameters of the foam [20]. These moduli can also be con¬ 
nected to the moduli of hber networks as has been done by [24, 25, 26]. But, for 
large compressions these moduli are no longer sufficient and we must use more 
sophisticated theories. 

Our goal in this paper is to further study the mechanical and structural re¬ 
sponses of fibrin to compression by a combination of experiment and theory. In 
particular, from precise dynamic rheological and microscopic measurements we 
have developed a model for the compression behavior of fibrin networks based 
on the theory of foams and other cellular solids. First, the compression stress- 
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strain curve was measured in experiment along with the loss and storage moduli 
at various strains with simultaneous 3D visualization of the fibrin network dur¬ 
ing the deformation. In addition to the non-linear behavior of storage and loss 
moduli, we were able to reveal the non-uniformity of the compressive deforma¬ 
tion with formation of a “compression front” or “phase boundary” along the 
axis of compressive strain (see Figure 1) [27, 28]. Next, we described how the 
non-linear rheological behavior could be explained by viewing the fibrin net¬ 
work as a “cellular” solid that could exist in two phases - the low-strain phase 
in which the fibers are mostly straight, and a high-strain phase in which the 
fibers are mostly buckled. In displacement-controlled experiments such a solid 
has a stress plateau when there is a mixture of both phases. The fraction of each 
phase at equilibrium is determined by a lever rule [29]. However, the storage 
and loss moduli at various degrees of compression in our fibrin networks were 
obtained in the experiments by performing small oscillations around a certain 
strain. These oscillations changed the fractions of the two phases but they were 
not necessarily quasi-static. We treated these two phases using a continuum the¬ 
ory of phase transitions in which a kinetic equation relates the rate of change 
of the fraction of the phases to the stress [30]. In this way we explained the 
storage and loss moduli of the networks while accounting for the contributions 
of the moving phase boundary as well as the viscoelasticity of each phase. 

In this paper we first provide the details of the experimental procedures and 
present a description of the obtained data on controlled fibrin compression fol¬ 
lowed by a comprehensive theoretical analysis. After a short summary of the 
well-known results for cellular solids we show how this approach can be applied 
to hber networks. We also describe the constitutive laws for hber networks 
under large compression that have been studied using experiments and theory 
for several decades [32, 33]. We finally show how these results can be gener¬ 
alized and combined with a continuum theory for phase transitions to obtain 
the storage and loss moduli of networks under compression. The novelties of 
this paper are (a) experimental visualization of a propagating phase boundary 
in compressed hbrin networks, (b) an explanation for the trends in storage and 
loss moduli of the compressed networks in terms of a theory of phase transitions, 
and (c) hrst application of foam theory to a bio-polymer network. Our analysis 
leads to some new experimentally testable predictions about the dependence 
of the mechanical properties of the clot on network parameters that have been 
summarized in the discussion section. 


2 EXPERIMENTAL PROCEDURES 

2.1 Fibrin Clot Formation 

Fibrin clots were formed by adding human a-thrombin (1 U/mL final con¬ 
centration) to pooled human citrated platelet-poor plasma in the presence of 
CaCb (40 mM final concentration). Before clotting, Alexa-Fluor 488-labeled hu¬ 
man hbrinogen (Molecular Probes, Grand Island, NY) was added to the plasma 
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Linear Phase ■ Densified Phase 



Figure 1: A fiber network under compression, (a) The top panel is a cartoon, 
schematically showing the high strain phase with a densified network on the 
right, the low strain phase with straight fibers on the left and a mixture of these 
two phases in the middle. Here w indicates the frequency of applied oscillation 
and s is the rate of change of phase fraction due to this oscillation. In the middle 
figure a sharp phase boundary separates the linear phase in the bottom from 
the densified phase on the top. The fraction of densified phase increases with 
increasing compressive strain, (b) The bottom panel quantitatively shows the 
tri-phasic stress-strain response of such a material assuming quasi-static loading. 
Theoretical curves are plotted according to Eqns. (3) and (8). We have used 
I = 1.32/im, d = 0.22/xm, v = 0.1^m“^, and the two coefficients n = 3, k = 0.5 
corresponding to a fibrin network. 
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samples (5% of total volume, 0.08 mg/mL final concentration) to visualize fibrin 
network structure in a fluorescence confocal microscope. The clots were allowed 
to form at room temperature for 90-100 minutes to ensure covalent cross-linking 
by factor Xllla. Ten fibrin clots from different plasma samples were prepared 
as described and used for compression studies with combined structural and 
rheological examination. 

2.2 Visualization of Fibrin Networks and Image Analysis 

The fluorescent confocal microscopy and strain-controlled oscillatory mea¬ 
surements were performed synchronously using a home-built “confocal rheo- 
scope” which combined Bohlin Gemini rheometer with Nikon Eclipse 200/Con- 
focal VT Eye microscope [18, 34]. The plasma clots were allowed to form be¬ 
tween horizontal rheometer plates separated by a distance of 650-750 /rm. The 
sample volume was 300-375 fA. The upper moving acrylic plate had a 20-mm 
diameter and the base motionless plate was made of a 22-mm-diameter round 
transparent microscope glass cover slip to perform confocal imaging of the net¬ 
work during rheological measurements. 

To improve the image quality, in a separate series of experiments uncom¬ 
pressed and compressed fibrin clots were analyzed using the Zeiss LSM510 
META NLO laser scanning confocal microscope to get higher resolution im¬ 
ages with z-stacks spanning the entire clot thickness. We used Plan Apo 40x 
water immersion objective lens (NA 1.2). A laser beam with wave length of 
488 nm was used for fluorescent imaging. The distance between slices in each 
z-stack was 0.5 ^m. The image resolution of each slice was 1024x1024 pixels. 

The samples were subsequently processed and analyzed quantitatively using 
the automated program as described in [36]. Briefly, a new two-stage method 
was used to extract the structure of the network in the 3D confocal microscopy 
images and to characterize the structure by vertices (end points and branch 
points of the network) and edges (fiber network segments) in a 3D graph. A 
detailed description of the algorithm can be found elsewhere [36]. 

2.3 Controlled Compression of Fibrin Clots 

For confocal microscopy there were two types of fibrin clots compressed using 
a computerized AR-G2 Rheometer (TA Instruments, New Gastle, DE). First, 
the 650-750-/im-thick fresh hydrated fibrin clots were formed directly between 
the rheometer plates and compressed vertically down to 1/10 of their initial 
thickness. 12 clots from various plasma samples were tested. The stepwise 
compression of the clots was performed in 25-/im and 50-/im steps at the rate of 
30 /im/s, as the downward force was applied by the upper rheometer plate on the 
top surface of the clot. In the other type of compression experiments the 150-/rm- 
thick fibrin clots were formed in a flow chamber built of a microscope glass slide 
(bottom) and a glass cover slip (top) separated by an adhesive and highly plastic 
putty, which permitted irreversible compression of the clot samples within the 
chamber placed between the rheometer plates. In both types of experiments the 
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fully polymerized and cross-linked plasma clots were compressed to a different 
degree followed by confocal microscopy. We typically waited 2 minutes after each 
step before performing confocal microscopy. A compressive strain (or degree of 
compression), e was defined as the absolute fractional decrease in fibrin clot 
thickness £ = \AL/Lq\, where AL = L — Lq, and Lq and L are the initial 
and reduced thickness dimensions of the uncompressed and compressed clots, 
respectively. 

2.4 Oscillatory Shear Rheometry of Compressed Fibrin 
Clots 

Changes of viscoelastic properties of fibrin clots during compression were 
measured by oscillating rheometry using a Bohlin Gemini rheometer (Malvern 
Instruments, Westborough, MA) with parallel plate geometry. Measurements 
were performed with a fixed oscillating 0.5% strain amplitude at a frequency 
of 1 Hz to produce a linear stress-strain response to imposed shear. The os¬ 
cillations were initiated immediately after each subsequent step of compression 
and lasted for 300 s. In each test, the storage (G') and loss (G") moduli were 
measured in combination with confocal microscopy of the fibrin clot structure. 
The elastic response of the deformed clot to applied shear was characterized by 
the storage modulus, G', which represents the stored energy and is defined as 
G' = (to/ 7 o) cos((5). Here, 70 is a (small) strain amplitude, tq is a shear stress 
amplitude, and <5 is a phase shift in the shear stress with respect to applied os¬ 
cillatory strain. The viscous response of the deformed clot to applied shear was 
measured by the shear loss modulus, G", corresponding to the energy dissipated 
as heat, and calculated as G" = (To/ 70 ) sin((I). 

3 EXPERIMENTAL RESULTS 

3.1 Compression-Induced Non-Uniformity of the Fibrin 
Network Structure 

Visual examination of the fluorescent confocal images of the fibrin clots dur¬ 
ing compression revealed marked non-uniformity of the network structure. Here 
we define the XYZ-coordinate system as Z-axis along the vertical direction that 
is opposed to the applied compression. The XZ-plane sections of the hbrin 
networks clearly showed the appearance and progression of a gradient in the 
fluorescence intensity in the direction of compression with a more or less sharp 
boundary between the areas with different intensities (Figure 2). To precisely 
quantify the non-uniformity of the network density, we computationally recon¬ 
structed and analyzed the 3D images of a 150-/im-thick uncompressed hbrin 
clot and the same clot after 20% and 50% compression by segregating each 
image into 14 equal horizontal sublayers spanning the entire height of the net¬ 
work. The thickness of the sublayers at various degrees of compression was 
different depending on the height of the clot. We quantihed the uniformity of 
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Figure 2: The “compression front” or “phase boundary” is formed in response 
to vertical deformation of a fibrin clot. To provide visual examination of the 
changing clot structure, the distribution of fluorescence intensity in the XZ plane 
is shown in the same fibrin clot, uncompressed (left), about 20% compressed 
(center), and about 50% compressed (right). As the degree of compression 
increases, the fluorescence intensity (reflecting the network density) displays a 
gradient along the direction of compression. Here the compressive strain e is 
a degree of compression defined as e = |AL/Lo|, where AL = L — Lq, and Lq 
and L are the initial and reduced thickness dimensions of the uncompressed and 
compressed clots, respectively. The dashed lines outline the top (I), middle (II) 
and bottom (III) layers of the compressed clot, respectively. A horizontal arrow 
indicates the position of the compression front. 


network densification as the node density in each layer and plotted it against 
the distance from the bottom of the clot compressed from above (Figure 3). 
Our results showed that during vertical compression of a fibrin network a much 
higher degree of node densification was observed closer to the top of a clot, 
while the node and network density decreased towards the bottom. The top, 
middle and bottom layers of the compressed clots are indicated in Figure 2 and 
Figure 3. As the degree of compression increased from e = 0 to e = 0.2 and 
e = 0.47, the node density of the upper network layers increased by a factor 
of 1.8 and 2.3, respectively. At the same degrees of compression (e = 0.2 and 
e = 0.47) the node density of the middle layers of the network increased by only 
1.1 and 1.2, respectively. The density of the bottom portion of the network did 
not change at e = 0.2, but increased by 1.2 for e = 0.47 as is apparent from the 
diffuse region near the bottom of the clot in Figure 2 (right panel). However, 
this does not contradict the presence of a moving compression front from the 
top of the clot. Since the width of the front given by the parameter c in Eqn. 
(1) is about 30/im and its center is located at Zq = 42^m, the bottom portions 
of the network likely experience a non-zero strain that pushes them against the 
glass surface causing a moderate increase in fluorescence intensity. 

To confirm the structural non-uniformities in fibrin clots produced by com- 
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pressive loads, we also complemented the 3D structural analysis of compressed 
fibrin networks with an experiment to follow the displacement of microscopic 
fluorescent beads embedded into the network in response to clot compression. 
2-/im polystyrene fluorescent beads were tracked during the compressive defor¬ 
mation in a clot volume of 35.8 x 35.8 x 25.5/im in the bottom portion of the 
clot. The compression was applied in 25/im steps. After each step we waited 2 
minutes before measuring the displacements of the beads. We found that the 
beads’ position changed non-linearly in the Z-direction as a function of compres¬ 
sive strain (Figure 4), indicating the non-uniform compression of the network. 
In the case of spatially uniform compression, one would expect to observe a lin¬ 
ear displacement of beads with compression irrespective of their initial position 
and the distance from the bottom of a clot. However, in our experiments the 
beads exhibited less than 20% relative displacements when the clot was exposed 
to compressive strains from e = 0 to £ = 0.6, but revealed a rapid descent toward 
the bottom of the clot at £ > 0.6. To be more specific, when the compressive 
strain reached 0.6, the absolute distances changed by 2,1.5 and 0.6y^m for the 
beads initially located at 14.7,11.1 and 2.7/rm over the bottom, respectively. 
At a compressive strain of 0.7 the bead, which was the closest to the bottom 
of the clot, stopped moving because it had reached the surface. As the strain 
increased, other beads moved downwards until they also reached the bottom of 
the clot. The trajectory of all the beads turned downward at a similar com¬ 
pressive strain because they are all located within 15/im of the bottom while 
compression is applied in 25^m steps and the width of the front is about 30/im 
(see Figure 2 and Figure 3). The nonlinear response of bead displacement to 
network compression indicated that the deformation of the network occurred 
non-uniformly with the top layers being compressed earlier and stronger than 
the lower network portions. As a result, the beads did not undergo large dis¬ 
placements at low strains but their vertical position started to change drastically 
at about e > 0.6 as the boundary between the relatively compact (upper) and 
loose (lower) portions of the network approached the bottom of the clot. 

Thus, both the structural analysis and bead tracking provide evidence for 
the formation and propagation of a “compression front” inside the network as 
more dense upper layers overtake the less dense bottom portions of the fibrin 
clot during progressive compression. Such a front can be represented by a strain 
profile given as: 

e = a——), (1) 

c 

z{Z) = d -I- (1 — a)Z — 6clog ^cosh (^?^——)^ > (2) 

where Z is the original height of a material point in the unstressed (reference) 
configuration, Zq is the location of the center of the front in the unstressed 
configuration, and z is the deformed height that is measured in the experiment. 
We got Eqn. (2) by integrating dzjdZ = 1—e with boundary condition z(0) = 0, 
as in our experiments. Capital Zq cannot be measured in our experiment but 
using Eqn. (2) the change in Zq can be obtained from the experimentally 
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Figure 3: (a)The node density of the fibrin network as a function of height 
from the bottom of the sample based on 3D reconstruction of the network, 
uncompressed (e = 0) and compressed to different degrees {e = 0.2 and e = 
0.47). The node density was normalized with respect to the node density of the 
uncompressed clot. Note that the node density increases from bottom to top. 
The curve is flatter near the top and bottom of the sample and has a larger 
gradient near the middle. This suggests that the fibers are buckled near the top 
and straight near the bottom. The symbols represent experimental data (MSD, 
n=4) and we have modeled the regime of high gradient as a compression front 
or phase boundary. The lines represent fitting curves made using Eqns. (1) and 
(2) together with the experimental data with fitting parameters a = b = 0.2, 
c = 30, and Zq = 42 for red curve, Zq = 100 for green curve respectively. Z^ 
represents the center of the phase boundary in the reference configuration. The 
decrease of Zq shows that the front moves downwards in response to increased 
compression. The top (I), middle (II) and bottom (III) layers of the compressed 
clot (e = 0.47) are separated by the vertical dashed lines that correspond to the 
layers shown in Fig. 2. Large circles with vertical arrows indicate the position 
of the compression front, (b) Experimental data of node density under large 
compression is fitted by Eqn. (7). This justifies our use of Eqn. (8) for the 
densified network. We have used I = 1.32/im, d = 0.22^m, v = 0.1/im“^. 
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measured changes in the location of the phase boundary as a function of 0 as 
explained below. The network density is proportional to 1/(1 — e) as discussed 
in section 4.1. Thus, we can plot the network or node density as a function 
of z. Since the minimum and maximum values of the strain in Eqn. (1) are 
a ± &, we expect that 2b should approximately equal the transformation strain 
or the difference in strains between the high-strain and low-strain phases at the 
plateau stress cto (see Figure 1). a should be the strain at the center of the 
front, c should be a measure of the width of the front, and d is obtained by 
enforcing z(0) = 0 (see analysis in Appendix). Keeping this in mind we chose 
reasonable values of a, b and c to ht the experimental data. The data and the 
fits are shown in Figure 3. We found that 2b « 0.4, as expected to be the 
transformation strain shown in Figure 1. Furthermore, Zq comes closer to the 
bottom of the sample when higher compression is applied, confirming that the 
“compression front” moves downward in response to the loading. 



Figure 4: The absolute vertical distance from the bead to the bottom of the 
clot along the direction of network compression for various locations in the 
same clot compressed from e = 0 — 0.9. The network is characterized by z - 
the bead distance from the bottom at a certain degree of compression and zq 
- the bead initial distance from the bottom of the network. 2/rm polystyrene 
fluorescent beads were tracked during the compressive deformation in a clot 
volume of 35.8 x 35.8 x 25.5/iTO in the bottom portion of the clot. Densihcation 
occurs at the bottom of the clot only at high compressive strain, likely coinciding 
with the arrival of the compression front at the bottom. Before the arrival of 
the compression front the density in the bottom of the clot remains almost 
independent of strain. This is what we expect in the plateau regime when the 
network consists of a mixture of two phases that accommodate the increasing 
compressive strain simply by advancing the compression front. 
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3.2 Non-linear Mechanical and Structual Response of Fib¬ 
rin Networks to Compression 

Structural alterations of fibrin network during compression underlie changes 
of shear viscoelasticity studied with rotational rheometry. The rheological mea¬ 
surements of fibrin clots with various initial shear storage and loss moduli were 
performed with two different compression steps (10 /rm and 25 ^m ) and at differ¬ 
ent degrees of compression. The experimental results are shown in Figures 5 and 
6 in combination with the theoretical curves. The normal stress, calculated as 
the force distributed over the area of the smaller rheometer plate was measured 
for each degree of compression, and presented as the plots of normal stress, a, 
versus compressive strain, e (Figure 5). The curves were almost linear at lower 
degrees of compression with a steep increase of the normal stress at higher de¬ 
grees of compression (e > 0.5). The evolution of the shear viscoelastic properties 
of fibrin networks in response to compressive deformation displayed a typical 
tri-phasic behavior (Figure 6). In the beginning of compression (e < 0.05), the 
shear elastic modulus of a fibrin network was almost constant followed by a 
long relative plateau up to about e = 0.7 — 0.8 compressive strain and finally 
dramatic increase at the highest compressive strains. The loss modulus followed 
the same trend (Figure 6). If the current rheological measurements are corre¬ 
lated with the structural analysis of the network dynamics performed earlier in 
[18], the results could be described as follows. First, a linear viscoelastic re¬ 
sponse to compression was observed, in which most fibers are straight. Then, a 
plateau regime follows, accompanied by an increasing number of buckled fibers 
while the stress remains nearly constant as strain increases. Finally, we see a 
regime in which network densification occurs with a stress-strain response that 
is markedly non-linear and dominated by bending of fibers after buckling and 
inter-fiber contact. In the next section we will quantify these observations using 
a model for foams. We emphasize that the experimental results reported here 
are substantially distinct from those reported earlier [18] because in addition to 
the confirmed non-linear response, we were able to reveal the non-uniformity of 
the compressive deformation with formation of a “compression front” or “phase 
boundary” along the axis of compressive strain (see Figure 3) [27, 28]. 

4 THEORETICAL ANALYSIS OF THE EX¬ 
PERIMENTAL DATA 

4.1 Short Review of Cellular Solids 

A foam is a porous low density structure consisting of cells or fibers made 
of a viscoelastic material. In our case these fibers are made of fibrin. If Eg is 
Young’s modulus of a single fibrin fiber, then the Young’s modulus and shear 
modulus of the network or foam consisting of these fibers is given by [20]: 

El = A.0^, (3) 
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where the important non-dimensional quantity </> = ttvcPI/A is the volume frac¬ 
tion of fibrin in the network, v, d, and I are respectively the number of hber 
segments per unit volume, fiber diameter and fiber segment length between 
branch points, (po is the initial fiber volume fraction before compression. These 
expressions are derived under the assumption of small strains so that the stress- 
strain relation of the foam can be approximated as linear. We will designate 
this phase of our fibrin network as the linear phase or the low-strain phase. 

When this network is loaded in compression fibers will buckle [20] when a 
critical force is reached. From Euler’s formula, the critical axial load that causes 
fiber buckling is 


p — 

J- r r — 


n^TT^EsI 

P ’ 


( 5 ) 


where Eg, I and I are the Young’s modulus, area moment of inertia, and the 
length of the fiber between branch points, respectively. The coefficient n is close 
to one and is determined by boundary conditions at the fiber ends. According 
to the analysis given in [20], n is taken to be 1 because the fiber lengh I varies 
over a range in random networks. We can use this expression for the critical 
load to estimate the stress (Tq at which buckling of hbers begins in the network. 
According to [20], since, (Jq oc Fcr/P oc Egl/l^, we could write it as 


(Jq — cEgcjyQ, 


( 6 ) 


where the coefficient c is assumed to be 0.04 in [20]. The strain corresponding 
to the critical stress in the low-strain phase is easily computed as ei = aojEi. 
Another way to estimate £i is to use Eqn. (5) - £i « Fcr/{AEs) where A is 
fiber cross-sectional area. Both these methods give an estimate of Si around 0.05 
and (To around 5 Pa as shown in Figure 1 for the parameters in the experiments 
shown in Figures 5 and 6. These ideas for estimating ctq can only be justified 
if the fibers in the network are not undulating due to thermal motion. The 
extent of thermal undulations depends on the persistence length of the fibers 
and their length between branch-points. The persistence length of hbrin varies 
over a wide range depending on the conditions (thrombin concentration, pH, 
calcilum concentration, etc.,[17]) under which a clot is made [21, 22]. In fact, 
[16] summarizes several reasons in addition to fibrin assembly conditions (such 
as, protein packing density, lateral binding between protofibrils, etc.) that 
result in the variation of Hbrin persistence length. For our fibrin networks, 
which have been derived from plasma clots, the average diameter of the fibers 
is about 220nm [18], which together with a Young’s modulus E = 5MPa [23] 
gives a persistence length EI/ksT at room temperature about Im. The length I 
between branch points is about 1 — 2/im, on average. Thus, thermal fluctuations 
of the fibers are negligible, consistent with the assumption for foams. 

When large compressive strain is applied to a fibrin network the total volume 
decreases considerably and the inter-fiber space is reduced. In this configuration 
a linear stress-strain relation is not applicable any more. The fibers are mostly 
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buckled and forced to touch each other. We call this kind of configuration the 
“densified phase” or the high-strain phase. This sort of densified network has 
been studied by Toll [33] who showed that the number of contact points W per 
unit volume is an increasing function of fiber volume fraction (j) and was given 

by 




16^ 

TT^ Cp 


( 7 ) 


where / is a scalar invariant of the fiber orientation distribution function and is 
equal to 7r/4 for three dimensional isotropic networks [33]. Following [20], in the 
densified phase we assume the fiber volume fraction (j) = (j)o / (1 — e). In other 
words, all the volume change is due to the reduction in the thickness of the 
network in the direction of the compressive force with no change in the cross- 
sectional area. Our experiments indicate that this is a good approximation. 
We compared Eqn. (7) to the experimental data on contact point density 
from [18] and found that it indeed varied quadratically with the fiber volume 
fraction (j) as shown in Figure 3. The stress-strain relation for such a network 
was proposed by van-Wyk [32] and Toll [33] in the following form: 


a = kE, (</)" - </)^), 


( 8 ) 


where the coefficient k is determined by the material and loading conditions. It 
is found to be less than 1, as in [35]. The exponent n was analytically derived 
to be 3 from Eqn. (7) for three-dimensional random networks in [33]. This has 
been confirmed by several experiments [35, 38, 39, 40]. Erom Eqn. (8), we can 
compute the local tangent moduli of the network in the densified regime. The 
tangent modulus at any strain in the densified regime is directly proportional to 
the storage modulus in the rotational rheometer oscillation experiments. Eol- 
lowing the linear theory for foams at small strains we assume that the ratio 
between this storage modulus and the local tangent modulus in the densified 
regime is 3/8 (see Eqns. (3) and (4)) and get the following expression for the 
shear modulus of the densified network at strain e:. 


G = 


^kE(j)l 

8(1-£)"' 


(9) 


Any contributions due to the Poisson effect in the densified regime are absorbed 
into the constant k that has been treated as a fitting parameter in the literature. 
This is also consistent with our experimentally motivated assumption above that 
vertical compression of the fibrin network does not change its cross-sectional 
area. For the plateau stress cr = ctq, we find that the Young’s and shear moduli 
are 


E2 


3kE(j)l 


( 10 ) 


G2 

where £2 is the strain at which ctq 


9kE(l)l 
8(1-£ 2 )"’ 
kEs{4>^-ct>l). 


( 11 ) 


13 



We will use a linearized version of the stress-strain relations (in the neigh¬ 
borhood of CT = CTo) in our subsequent analysis. These are summarized for the 
low- and high-strain phases as follows. 


f a I El cr < (To (e < £i) 

\ ujEiEeT cr > (To (e > £2) 


( 12 ) 


where cto is the critical stress for buckling, as discussed previously. At this 
stress the straight and densified phases co-exist. Ei,E 2 are Young’s moduli in 
the straight and densified phases respectively from Eqns (3) and (10). Et is the 
transformation strain which can be estimated by constructing the tangent to 
the stress-strain curve at £ = £2 and reading off the intercept on the e axis as 
shown in Figure 1. If we assume that the process of loading is quasi-static so 
that the sample stays in equilibrium during the compression process, then the 
compressive strain e in the plateau corresponding to the stress (Tq is obtained 
simply by changing the fraction of fibers buckled. If we denote the fraction of 
fibrin in the high strain phase as s € [0,1], then total strain can be estimated 
from the insight that the sample consists of two ‘springs’ in series consisting of 
the straight and densified network respectively as observed in Figure 2. 


£ — (1 — s) — -I- s -I- £t^ . (13) 

This completes the description of the stress-strain curve for a fibrin network 
under compression. As an example we have plotted this stress-strain curve for 
a fibrin network with network parameters experimentally measured in Kim et 
al. [18] in Figure 3. We have used values of the fiber Young’s modulus Eg 
= 5 MPa [9, 10, 41]. The stresses in Figure 3 are of the same magnitude as 
the measured storage moduli in Figure 6. However, the stress strain curve in 
Figure 5 does not have an initial linear regime for small strains even though 
it seems to have a flat plateau regime for intermediate strains and a regime 
with steeply increasing stress for high strains. Furthermore, the shear moduli 
measured in Figure 6 are at least three orders of magnitude lower than those 
seen in the normal stress data plotted in Figure 5. This indicates that Figure 5 
is not a true reflection of the stress strain response of the network; rather it is 
the result of water squeezing out of the network in response to compression. To 
understand why this is the case we must account for the strain-rate dependence 
of the mechanical response of these fibrin gels. 


4.2 Strain Rate Dependence 

The fibrin networks that have been tested in our experiments really are 
gels that contain a large amount of water. As the gel is compressed, water is 
squeezed out. However, as it is sheared, water is not squeezed out because the 
total volume does not change. Hence, it is expected that the measured normal 
stresses in compression will depend on the applied strain rates. This effect is 
well-known as poro-viscoelasticity in foams and other cellular solids [20]. An 
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Table 1: Fitting parameters for each experiment. 


Individual group 

Sample 
height (mm) 

strain 

rate(s“^) 

Fiber 

length(/im) 

Black circle 

550 

0.055 

0.5 

Magenta triangle 

650 

0.046 

0.3 

Black diamond 

650 

0.046 

0.27 

Blue circle 

750 

0.04 

0.5 

Green diamond 

650 

0.046 

0.4 

Red square 

750 

0.04 

0.3 


Individual group 

4’o 

k 

£20 

M (Pa”^s-i) 

WTo „(10 ^®J-s) 

Black circle 

0.0045 

0.025 

0.82 

0.004 

2 

Magenta triangle 

0.0045 

0.0025 

0.8 

0.005 

1.2 

Black diamond 

0.0035 

0.0025 

0.85 

0.008 

1.7 

Blue circle 

0.0035 

0.333 

0.7 

0.008 

4 

Green diamond 

0.004 

0.333 

0.6 

0.008 

6.5 

Red square 


Densification not observerd 


analytic relation for the stress as a function of the applied strain rate is given 
as [20] 


Cue f D\ 
= —e\j) 


(14) 


where /x is dynamic viscosity of the fluid, e and e are the compressive strain and 
its rate, D is horizontal dimension of the foam sample, and I is the cell edge- 
length of the foam. The coefficient C is about unity. As shown in Figure 5, 
the measured normal stress during compression can be captured by Eqn. (14) 
with the parameters /i = 0.001 Pa-s (for water), D = 22 mm (known from our 
experimental set-up) and I as the only htting parameter for each experiment. 
The strain rate for each experiment is different because the height of the sam¬ 
ples varies while the rate of compression (30/xm/s, see section 2.3) remains fixed 
for each experiment. The fitted values of I appear in Table 1. They are within 
the range of variation of fiber lengths between branch points seen in experi¬ 
ments [18]. We see that most of the normal stress measured in the experiments 
is due to the expulsion of water from the fibrin gel. Hence, the storage and loss 
moduli measured in our rheometer experiments and shown in Figure 6 cannot 
be extracted from the stress-strain curves shown in Figure 5. 

From Eqn (14), we can approximate the viscosity of the gel for small strains 
(infinitesimal values of e) as 

’7 0c/x(y) (15) 


However, the viscous losses due to this term cannot be measured in the rotational 
rheometer experiments because the derivation of this expression assumes that 
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the volume of the network is changing while the oscillatory experiments are 
performed in shear which involves no change in volume. Therefore, the loss 
modulus in the low strain phase should only come from the liquid viscosity 
itself. For water, this is merely /x = 0.001 Pa-s, which is much smaller than 
the value obtained from (15). Therefore, we treat the fibrin network in the low 
strain phase as a purely elastic solid with constant modulus. 

Another important network property that is related to the strain rate is 
the inter-fiber friction. This is particularly important in the densified phase 
in which a large number of contacts are created between fibers. Fibers sliding 
against each other while in contact will cause dissipation in the network. In 
order to estimate the viscosity associated with this process we appeal to an 
elegant calculation in [43]. This book gives an explanation for the molecular 
basis of “viscosity” due to the forming and breaking of bonds as molecules slide 
past each other. The expression is: 

V = wNcTon, (16) 

where w is the inter-molecular bond energy and Ton the average life-time of a 
bond. An expression for No as a function of fibrin volume fraction is already 
given in Eqn. (7). However, w and Ton are difficult to determine for fibrin fibers. 
Some work along this direction has been done for estimating the viscous losses 
in microtubule networks [42], but here we treat the product WTon as a fitting 
parameter. In summary, we model the densified network at ct = (Tq (for small 
strains) as a Kelvin-Voigt solid [44] in which the stress-strain relation is given 
as: 

a = E2{£ - £t) +'qe, (17) 

where the parameters E 2 and q are taken from Eqns (10) and (16). In the 
rotational rheometer experiments the loss modulus is related to this viscosity 
as: 

G" = quj, (18) 

where oj is frequency of oscillation. The loss modulus for both the low- and high- 
strain phases can be determined using this formula if we use q = q = O.OOlPa-s 
(of water) for the low-strain phase and Eqn. (16) for the high-strain phase. In 
order to estimate the loss modulus in the plateau regime of the stress-strain 
curve we need to account for the rate at which the fractions of the two phases 
evolve. This is explained in the next section. 

4.3 Phase Boundary Mobility 

In the oscillation experiments, the fraction of each phase depends on both 
applied strain and strain rate. As shown in Eigure I, the phase boundary is 
expected to oscillate as the applied loading oscillates. We treat this process 
in a one-dimensional model as follows. After pre-compressing the sample, we 
assume that the oscillation is also in the axial direction. This is different from the 
experiments in which an oscillating shear strain is applied on the pre-compressed 
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networks. However, since the oscillating shear strains applied in the experiments 
are small, we can assume that they will give storage and loss moduli that are 
proportional to those obtained in our calculations. 

During the oscillation process the sample may not be in equilibrium and thus 
the stress may not be (Tq. There will be oscillation in stress, phase fractions, as 
well as strain. We assume that the fraction s is determined by the kinetics of 
transformation as in [30]. This means that the rate of change of each fraction 
is given as a function of stress cr. To simplify the problem we consider the 
particular kinetic equation [31]: 

s = $ (cr) = M {a — (To), (19) 


or, 

M{a — ao)dt = M / {a — ao)dt, (20) 

Jo 

where M is a time-independent mobility parameter determined by the network 
that could be regarded as a fitting parameter in our theory and sq is the fraction 
of high strain phase before oscillation is applied. In the experiment we apply a 
harmonic oscillating strain to the pre-compressed network as: 



e = eo -I- a sin wt. 


( 21 ) 


where Eq is an initial strain in the plateau regime of phase transition before 
oscillation is applied. Taking account of the strain rate dependence we rewrite 
Eqn. (13) as: 


£ = (1 - s)£l -I- S£2, 

a = EiSi = i?2(£2 — £t) + VS2- 


( 22 ) 


Now, recall that the expression for equilibrium stress (Jq and strain Sq are given 
by: 


£o — (1 — so)£io + •So£207 
Co = EiSio = i?2(£20 — £t)- 


(23) 


where the subscript 0 denotes the quantity in initial compression equilibrium 
before oscillation. Subtracting Eqns. (23) from (22) we get: 


asincot = (1 - so)a; + soy + (s - so )(£2 - ei), (24) 

CT - (To = Eix = E 2 y + Tjy, (25) 

where x = £i —eio and y = £ 2—£20 are the small oscillating components of strain 
in each phase. Since the oscillation in strain is small, we make an approximation 
in Eqn. (24) by replacing £2 — £1 with £20 — £io- Then, plugging Eqn. (20) and 
(25) in (24) we get: 


a sinwt = (1 — sq) 


E2y + rjy 
E^ 


Soy 


-kM(£2o-£io) [ {E2y + vy)dt. 
Jo 


(26) 
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Carrying out the integral of y and rearranging terms, we find that 


asinujt = Ay + By + C / ydt, (27) 

Jo 

where the three dimensionless quantities A, B, and C are given by: 

E 2 

^ = (1 — ^0 )^—f So + Mri{e 2 o ~ £10)) 

C = ME2{e2o — £io)- 


Let us assume that this equation has the solution: 


y = U sinojt+ V cosojt, (29) 

where U and V are at present unknown. Plugging this into Eqn. (27) and 
solving for the two coefficients U and V gives: 


^ “7l2 + (Bw-C/w)2’ 

Bu) — C ju! 

^ ^ ““^2 + (Bw- (7/^)2 ■ 

Plugging this back into Eqn. (22) we get the expression for stress as: 


(30) 


cr = (To + £^ 2 y + yy 

= (Jo + a{E 2 U — rjVuj) sinwt + a{E 2 V + rjUuj) coscot. 

Therefore, the storage and loss modulus in the plateau regime (i.e. a = uq) are 
respectively 


G' = E 2 U -yVoj 
G" = E2V + yUio 


E 2 A + rjBuj'^ — T]C 
A^ + {Bi^-G/ujf' 

—E 2 {Buj — C/w) + yAu) 
A^ + (Buj - Cjijuf 


(32) 


In order to confirm that the above framework predicts the correct storage 
and loss moduli, we compared quantitatively our theory with the experimental 
data. The input parameters of the network in the above computations were 
taken from the experiments shown in Figure 6. Combining Eqns. (4), (9) and 
(32) together, we obtain the theoretical curve in each individual experiment in 
Figure 6. We take the value of branch point density and fiber density exactly 
as measured in experiment as 0.04/xm“^ and 0Ay,m~^ respectively. The other 
fitting parameters for each experiment are reported in Table 1. The trends in 
storage and loss moduli are captured quite well with this model. Note that the 
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experimentally measured average fiber length and average fiber diameter are 
1.32 fj,m and 220 nm respectively [18], which result in a fiber volume fraction 
of <j)o = TTvd^l/A = 0.005. The (po values obtained by our fits are slightly below 
this value but well within the range of variation of (pQ in the experiments [18]. 
Also note that the value of constant k in the densifed phase varies inversely with 
the strain £20 where the phase transition finished. We have not yet been able 
to rationalize this. 

We will now connect our results for the variation in moduli to Figure 1 which 
schematically shows three distinct regimes in the stress-strain curve. To this end 
we have replotted all the data for storage and loss moduli in Figure 7 against 
the logarithm of the strain e. The red and blue curves in the upper and lower 
panels are fits to the average of all the storage and loss moduli, respectively. 
The strains £1 and £2 have been indicated by dashed lines to demarcate the 
three regimes. In the initial linear regime, £ < £ 1 , the moduli G' and G” 
are constant. In the plateau regime, £1 < £ < £2 the network accommodates 
more strain by changing the fractions of the straight and densified phases of 
the networks through the motion of the phase boundary. G' decreases with 
increasing strain in the plateau regime because the slope of the stress-strain 
curve of the densified phase at £ = £2 is smaller than the slope for the straight 
phase at £ = £ 1 . For £ > £2 densification occurs and both moduli increase 
steeply. 

5 DISCUSSION 

We have shown in this paper that fibrin networks are natural cellular mate¬ 
rials. The theoretical description of the fibrin mechanical behavior is novel and 
treats fibrin networks in the context of a broad class of natural and synthetic 
foam-like materials which share some key features. Among these is the tri-phasic 
stress-strain response with a rarefied low-strain phase, a densified high-strain 
phase, and a plateau phase consisting of a mixture of these two. Like most 
foams we also see a moving phase boundary separating the two phases when 
fibrin networks are compressed [27, 28]. The phase boundary and several other 
features quantifying the structural non-uniformity of fibrin networks under com¬ 
pression were revealed by our 3D microscopy studies. Because of the complexity 
of the analysis, non-uniform fibrin networks under loads have not been previ¬ 
ously quantified. However, it is clear that under physiological conditions fibrin 
networks are frequently exposed to dynamic forces, which would result in non¬ 
isotropy and structural heterogeneities. Since in many cases, fibrin networks 
participate in the permeation of blood clotting factors [45], variations in fib¬ 
rin matrix density can affect the spatial-temporal distribution of factors and 
hence, alter their biochemical interaction. Our results indicate that even a uni¬ 
form compression can produce structural gradients in fibrin networks, therefore, 
suggesting a non-uniform alteration of fibrin matrix transport properties with 
compression. 

It is also worth noting that the mechanical response of fibrin networks un- 
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der compression is dominated by fiber buckling and subsequent bending and is 
therefore non-afhne. This is in contrast to the stress-strain response in tension 
and shear where the assumption of affine deformations worked quite well [9, 10] . 
This assumption is particularly good for rubbers and other polymeric materials 
[11, 12], and it works for fibrin networks under tension because the unfolding 
of the proteins results in cross-linked polymer chains much like those in rubber. 
It works well for shear of a network of thin fibrin fibers because such a net¬ 
work is not very different from cross-linked fluctuating polymer chains [14, 16]. 
However, in compression the initial linear regime is governed by fiber bending 
instead of fiber stretching. Thus, we can no longer assume affine deformations. 
This is apparent in the non-linear dependence of the initial Young’s modulus 
in compression Ei on the fibrin volume fraction (j) (see Eqn. (3)). The depen¬ 
dence would be linear if the deformations were affine as in tension and shear. 
Since our model predicts that Ei oc where v is fiber den¬ 

sity, d is average fiber diameter and I is average fiber length. Our prediction 
can be tested in experiments since v, d and I can be controlled by changing 
the polymerization conditions of the fibrin networks. In the densified network 
buckled fibers are in contact, so the assumption of affine deformations is not 
expected to work. Indeed, the stress-strain response depends on in a complex 
way that has been quantified both with experiments and theory in the textile 
literature [32, 33]. We have shown in this paper that the same ideas can quan¬ 
titatively capture the compression response of densified fibrin networks with an 
appropriate choice of fiber geometrical and material parameters. For example, 
at large compressive strains our model predicts that the tangent modulus in the 
densified phase E 2 oc <Pq/{ 1 — e)'^ where (pQ is the initial fiber volume fraction. 
This prediction can be tested since v, d, I and e can all be controlled in ex¬ 
periments. Similarly, our physically motivated constitutive law Eqns. (16) and 
(7) that accounts for the inter-fiber friction and captures the viscous response 
of densified fibrin networks can also be tested in experiments. We have already 
demonstrated the importance of the strain-rate dependence in the compressive 
response of fibrin networks by using poroelastic constitutive laws that govern 
the flow of water in cellular solids. 

In this work the structural alterations in response to compression were stud¬ 
ied at the microscopic level of a network. We have shown that structural non¬ 
uniformity manifests itself as a “compression front” or “phase boundary” sep¬ 
arating regions in the network where fibers are straight from another regime 
where the fibers are buckled. While such a front has been visualized in many 
macroscopic [27, 28] and microscopic foams [19] we are not aware of attempts 
where its motion is treated using a theory of phase transitions [31] to obtain the 
storage and loss moduli. We have demonstrated in this work that a continuum 
theory of phase transitions can predict the rheological properties of fibrin net¬ 
works under large deformations if we properly account for the dissipative motion 
of the phase boundary. It is likely that a similar model can be applied to other 
foams or foam-like materials since propagating interfaces have been visualized 
in them. The idea of a phase transition was also used in our earlier work on the 
tensile behavior of fibrin networks [9, 10]. The difference here is that we have 
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no evidence of forced molecular unfolding of the protein. However, fibrin has 
a complex multidimensional structure and it has been shown that compression 
of a hydrated fibrin clot is accompanied by changes in the a-helical coiled-coils, 
namely the transition of a-helix to /?-sheet revealed using FTIR spectroscopy 
[46]. However, a comparative analysis shows that the changes in the secondary 
structure of fibrin were observed at compressive stresses (> 25MPa) that are 
orders of magnitude higher than those applied in our work (Figure 5). It is 
likely that molecular structural transitions observed in a protein at very high 
degrees of compression results from its mechanical destruction, not reached in 
our present experiments. The aggregate of data suggests that at the degrees of 
deformation, not leading to complete disorganization of a filamentous network, 
compressive load is accommodated mainly by the structural changes in fiber 
arrangements at the network level. These changes include buckling, stretching, 
and formation of oblique contacts of individual fibers as well as overall network 
densification. It is noteworthy that the structural changes underlying compres¬ 
sive deformation of the networks are spatially non-uniform and form a stepped 
gradient along the axis of compression. Although we observed a single compres¬ 
sive front in our experiments that originated at the boundary of our sample, 
this is by no means necessary. Fronts can be nucleated in the interior of the 
network at a site of stress concentration or low network density. This prediction 
can also be tested since it is possible to visualize moving fronts as we show in 
our experiments. 

Fibrin deformability has become a rapidly developing area not only in the 
study of the mechanics of filamentous networks, but also in various practical 
biomedical applications. First, studies of compression of fibrin networks provide 
a basis for understanding what happens to salutary hemostatic clots and patho¬ 
logical obstructive thrombi in vivo, when they are subject to natural deforma¬ 
tions. There is a number of conceivable patho-physiological conditions, in which 
fibrin networks can experience compressive loads: vasoconstriction, arterial hy¬ 
drodynamic blood flow, clot contraction, pulsation of a vessel wall or aneurism, 
muscular pressure, cardiac beats for intracardiac thrombi, lung motion during 
breathing, etc. Next, fibrin has been widely used as a sealant in surgery and its 
hemostatic efficiency strongly depends on the ability to withstand mechanical 
loads, including compression. Finally, fibrin has been well-known for various 
tissue engineering applications in combination with cells and with other mate¬ 
rials to replace damaged or diseased organs and tissues [47]. Mechanical and 
structural insights into fibrin deformability in combination with a theoretical 
framework of studying fibrin mechanics provide a foundation for developing 
new treatment modalities. In particular, modeling of fibrin compressibility us¬ 
ing the theory of foam mechanics opens an avenue for physically treating fibrin 
clots and thrombi as well as for developing new fibrin-based biomaterials with 
predictable properties. 
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Appendix 

Structure of the Compression Front 

In our analysis above we assumed that the compression front is sharp. How¬ 
ever, our experiments reveal that it has a thickness over which there is a large 
gradient in the strain proHle. Fortunately, these descriptions are connected and 
have been discussed at length in the literature [30]. Recall that we Htted the 
strain profile in our compressed sample using Eqn. (1). In the following, we 
briefly describe why we use the strain profile given by Eqn. (1). We follow 
the analysis in [30]. Our network can exist in multiple phases at the same 
stress. Such materials have an up-down-up type stress-strain curve that can be 
described using 

cr - CTo =-a(£ - £*)-I-/3(e - e*)^, (33) 

where a is the stress, s is the strain, and a, /3, ctq and e* are material parameters. 
If we remember that cr = , where W (e) is a stored energy function, then it 

is not difficult to see that IF is a quartic function of e with two wells. The well 
at low strains corresponds to a network in which the fibers are straight, while 
the well at high strains corresponds to a densified network with buckled fibers. 
Without loss of generality we will now refer to cr — (Tq as cr and to e — e* as e. 
As an example, cro = 5 Pa, e* = 0.32, a = 51 Pa and j3 = 651 Pa captures the 
compression curve for Eg = 5 MPa in Figure 3. In Figure 8 we have plotted the 
Gibbs free energy density = W — as) landscape using the above parameters 
for three different values of cr to show which phase has lower energy at each 
stress. Since the motion of the compression front causes energy dissipation 
we augment this constitutive equation with a linear viscosity term and since 
interfaces cost energy we add a linear strain gradient term as follows [30]: 

a = ae-j3e^ + pv--pX — , (34) 

where p is the density of the material, v is a, viscosity and A accounts for the 
energy required to create an interface. We assume that the compression front 
travels at velocity v and define a new variable ^ = Z — vt. Since the stress cr in 
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our sample is constant we wish to find a profile for s(z, t) that satisfies 

s ds 

+ pvv—+ ae - + a = Q, (35) 

a4 

with remote conditions e —>■ e+ as ^ oo and £—>■£“ as ^ — oo. It is easy 

to show that 

£(z,t) = a + 5tanh(-- — -(36) 

c 

gives such a profile for a particular choice of a, b, c and v that depend nonlinearly 
on a, P, V, A, p and cr. In particular, o is a solution to the cubic CT + 2aa —8/3o^ = 

0, 6 = ^j — ia?, c = ^ and v = We recover Eqn. (1) by setting 

V = 0, justifying its use to fit the strain profile in our compressed samples. Note 
that in accordance with our assumption in Eqn. (19) the front velocity v in 
this augmented theory depends on ct. A sharp front corresponds to the limit as 
c —>■ 0. In this limit the dissipation caused by the motion of the front can be 
captured by a kinetic equation like Eqn. (19) as in [30]. 
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Addendum 


After the publication of our article [1] online it was brought to our attention 
that there is some overlap in its content with an article published by Lakes et 
a/. [2]. Unfortunately, we were unaware of this paper and did not cite it in our 
article. After carefully reading their paper, we provide a summary of the points 
of similarity and differences between the two papers below. 

1. Both papers study, by a combination of experiment and theory, the com¬ 
pression behavior of fibrous materials. They both show that the measured 
stress-strain curve has three regimes - a linear regime, plateau regime, 
and densified regime - typical of foams [3]. However, the Scott industrial 
foam studied in [2] is a macroscopic reticulated open cell synthetic foam 
with a dehnite pore size, while the fibrin studied in [1] is a microscopic 
isotropic network of natural fibers with a distribution of lengths and pore 
sizes that has not been examined in the context of theories for foams and 
cellular solids. 

2. Both papers show experimentally that in the plateau regime regions of 
densihed material (with inter-fiber contact) co-exist with rarefied material, 
and the fraction of densified material increases with increasing strain in 
the plateau regime. In [2] the bands of localized compression observed in 
experiments were viewed as regions of discontinuous strains alternating 
between two distinct values A+ and A_. In [1] we visualized a moving 
“compression front” with densified material on top and rarefied material 
in the bottom. In addition, we showed that this front starts at the top 
of the sample and reaches the bottom at compressive strains greater than 
0 . 6 . 

3. Both papers argue that the plateau in the stress-strain curve can be the¬ 
oretically explained by accounting for the co-existence of the rarehed and 
densihed phases of material at the same stress a, but distinct values of the 
stretch A in each phase. Lakes et al.[2] argued that their strain disconti¬ 
nuities are analogous to phase boundaries in continuum theories of phase 
transitions, and that minimization of total potential energy in the pres¬ 
ence of these phase boudaries leads to the stress plateau. Without prior 
knowledge of the work in [2] we explained the plateau in the stress-strain 
curve of compressed hbrin by the quasistatic motion of the “compression 
front”, or phase boundary, visualized in our experiments. 

4. Both papers point out that some of the energy dissipation observed in 
experiments can be partly attributed to the dissipative motion of phase 
boundaries. While this was a qualitative statement in [2] we have quanti¬ 
tatively deduced a mobility parameter for phase boundaries in compressed 
fibrin networks by fitting our measurements of the storage and loss mod¬ 
uli at various compressions with analytical expressions derived using a 
continuum theory of phase transitions [4]. 
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5. Both papers point out that the phase boundary is not a sharp interface 
but has a thickness over which there is a large gradient in the strain [1]. In 

[2] it was observed that the width of the bands depended on cell size. The 
authors mentioned that the length scale of the bands could be accounted 
for via models which incorporate strain gradients. In [1] we quantitatively 
fit the strain profile observed in our experiments with analytical expres¬ 
sions derived from an augmented theory of phase transitions incorporating 
strain gradients [4]. We show that the thickness of the phase boundary 
is many times larger than the average length of a fiber between branch 
points. 

The results of [2] are important and relevant to the experimental and theoret¬ 
ical study of foams and our paper, and we duly cite the article through this 
addendum. 
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Figure 5: Fitting of the experimental data of normal stress according to Eqn. 
(14). Different symbols and colors are experimental data corresponding to dif¬ 
ferent clot samples of varying height as shown in Table 1. Since the rate at 
which they are compressed is 30/im/s for all samples the strain rate for each 
of them is different. Lines show model fits with parameters as given in Table 
1. Most of the normal stress measured in the compression experiments on our 
fibrin gels is due to the escape of water in response to compression. For this 
reason we cannot extract the storage and loss modulus of the networks from 
these stress-strain curves. 


29 




















C5 


10 


10 




j I 5 j f 


1 5 


0.2 0.4 0.6 0.8 1 

Compressive Strain e 


Figure 6: Fitting of the storage moduli and loss moduli according to Eqns. (4), 
(9), (32), (16), (18), and (32). We have used fitting parameters for each individ¬ 
ual experimental group as given in Table 1. Note that there is no densification 
regime in the last group. One possible reason could be the failure or damage 
of the network under large compression. However, the data for strains smaller 
than 0.8 in the last group are consistent with the other data sets. 
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Figure 7: Experimental data of storage moduli and loss moduli versus logarithm 
of strain, showing three different regimes of fibrin network mechanics as indi¬ 
cated schematically in Figure 1. In the linear regime(0 < e < 0.05), the moduli 
are constant; in the phase transition regime(0.05 < e < 0.7), they are decreasing 
due to progressive fiber buckling; in the desification regime (0.7 < £ < 1), they 
are steeply increasing due to bending of the buckled fibers and increased inter¬ 
fiber contact. The theoretical curves passing through the experimental data 
are plotted using Eqns. (4), (9), (32), (16), (18), and (32) with (po = 0.0045, 
k = 0.025, £20 = 0.7, M = 0.004Pa"^s-i, ujTon = 2 x 10-i®J • s. 
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Figure 8: Multi-well Gibbs free energy landscape for various stresses. At the 
plateau stress a = 5Pa both wells are at equal height, meaning that the straight 
and densified phases of the network can co-exist. For lower a the straight phase 
has lower energy and for higher a the densified phase has lower energy. 
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